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Oxygen isotope effect on structural parameters and on the thermal motion in La 2 Cu04 
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We studied the oxygen isotope effect (OIE) on the zero point motion and the thermal motion 
as well as on structural parameters in La2Cu04 by means of high-resolution neutron diffraction 
experiments. We found a negative OIE on the lattice parameters (-0.01%). The OIE on the isotropic 
thermal parameters turned out to be positive for lanthanum, copper and negative for the oxygen 
atoms, respectively. The Rietvield refinement of the anisotropic thermal parameters of the different 
directions yielded the same sign of the OIE for each ion. Our analysis shows that the OIE on 
isotropic thermal parameters is most pronounced for La at 15 K (up to 100%) which we found to 
be originating mainly from the thermal motion in z-direction as determined from the refinement of 
the anisotropic thermal parameters. 

PACS numbers: 74.70.Dd, 74.62. Fj, 74.25.Ha, 83.80.Fg 



I. INTRODUCTION 

Possibly the largest enigma in contemporary solid state 
physics is how doping holes into the antifcrromagnetic 
planes can lead to high transition temperature supercon- 
ductivity in the cuprates. Early after their discovery, 
Anderson proposed an explanation in terms of resonat- 
ing valence bonds (RVB)i based on the fact that the 
parent compound exhibits long-range antiferromagnetic 
order—. Indeed quantum effects seem to play an essen- 
tial role in the cuprates: For instance the observed zone 
boundary spin- wave dispersion in L^CuO^ may be a 
fingerprint of an intrinsic quantum effect typical for a 
2D quantum Heisenberg antiferromagneti The impor- 
tance of quantum fluctuations should also be reflected in 
zero-point motion effects. In the case of Ceo the extent 
of zero-point motion effects turned out to be significant 5 
in a sense that they may be responsible for a substan- 
tial renormalization of electron-phonon coupling. In the 
lanthanum cuprates, this point has not been addressed 
experimentally so far. In particular no information about 
the oxygen isotope effect (OIE) on the zero point motion 
is available. 

Isotope effects on electronic properties in the cuprates are 
well known, the superconducting transition temperature 
decreases with oxygen isotope substitution^ whereas the 
pseudogap temperature raises as has been demonstrated 
e.g. by means of neutron cyrstal-field spectroscopy 7 . Re- 
garding isotope effects on magnetic properties only little 
is known so far. Recently Khasanov et al.^ performed 
a detailed muon-spin rotation and magnetization study 
of the isotope dependence of magnetic quantities, they 
found that the antiferromagnetic ordering and spin-glass 
ordering temperature in Yi_ 2; Pr 2 ;Ba2Cu307_ 1 5 exhibit a 
large OIE in the regime where superconductivity and an- 



tiferromagnetic order coexist. A huge OIE on the spin 
glass temperature was also found in Mn-doped LSCO at 
low doping while the isotope effect is small in Mn-frcc 
samples 9 . These unusual effects could arise from the iso- 
tope dependent mobility of the charge carriers 1 ^. As for 
the undoped compound La2Cu04, the Neel temperature 
Tjv slightly decreases upon oxygen isotope substitution^ 1 -. 
These findings were assumed to originate from structural 
changes upon oxygen isotope substitution — . But so far 
only the isotope effect on the lattice constants and the fol- 
lowing orthorhombicity have been measured. Thus there 
is need for a careful neutron diffraction study to extract 
the OIE on structural parameters and on thermal motion 
which might be relevant for J, the coupling constant of 
the in-plane antiferromagnetic superexchange interaction 
between nearest- neighbor Cu 2+ spins. 
To do so we perforemd a high-resolution neutron diffrac- 
tion study. Neutron diffraction is an excellent tool to 
study thermal and zero point motion since the neutron 
cross section is directly proportional to the Debye- Waller 
factor. 



II. EXPERIMENTS 

Polycrystalline samples of La2Cu04 were prepared us- 
ing conventional solid state reaction. Oxygen isotope 
exchange was performed using the procedure described 
previously^. The oxygen content was determined by 
thermogravimetric hydrogen reduction! 4 . Magnetization 
measurements were performed using a Quantum De- 
sign MPMS in fields ranging from 20 mT to 6 T at 
temperatures between 4 and 300 K on La2Cu 16 04 and 
La 2 Cu 18 4 . 

Subsequent neutron diffraction experiments were per- 
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FIG. 1: Susceptibility vs. temperature for La2Cu 16 04 (solid 
circles) and La2Cu 18 04 (open circles). The inset shows the 
derivative of the susceptibility. 



formed on the high-resolution diffractometer HRPT 15 at 
SINQ 16 located at the Paul Scherrer Institute, Switzer- 
land. The experiments were carried out at a wavelength 
A = 1.1545 A. The La 2 Cu 16 4 and La 2 Cu 18 4 sample 
were each contained in a 8mm-diameter Vanadium con- 
tainer which was then mounted into a closed cycle re- 
frigerator in order to reach temperatures between 15 K 
and 290 K. High statistics data were taken at 15 K and 
290 K (3.5 • 10 7 counts), whereas points at temperatures 
15 < T < 290 K were obtained at intermediate statistics 
(2.5 • 10 6 counts). 
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FIG. 2: Neutron diffraction data at 15 K for La2Cu O4 (up- 
per panel) and La2Cu 18 04 (lower panel). The solid black 
line corresponds to a Rietvield refinement as explained in the 
text. The difference between measured and calculated profile 
is shown below. Tick marks at the bottom represent the po- 
sition of allowed Bragg reflections in the Bmba space group. 



III. RESULTS 

A. Macroscopic measurements 

Susceptibility x{T) versus temperature is shown in Fig- 
ure Q] in an applied field of 0.1 T for the 16 O and the 
18 sample. A clear peak in the susceptibility is seen at 
261.4 ±0.1 K and 263.8 ±0.1 K for the 16 and the 18 
sample, respectively. The derivative dx/dT is shown in 
the inset of Fig. Q] 



B. Neutron diffraction 

The intensity pattern from neutron diffraction were re- 
fined using the orthorhombic nonstandard space group 
Bmba which is isomorphic to the standard space group 
Cmca (No. 64.), with atomic positions 01= (1/4, 1/4, z), 
02= (0,y,z), La= (0,y,z) and Cu= (0,0,0). The oc- 
cupancy of all atoms was set to the full value, the oxy- 
gen content was fixed to 4.0. Isotropic and anisotropic 



Debye- Waller temperature factors (short: thermal pa- 
rameters) were refined for each data set. The results for 
the lattice constants, atomic positions and Cu-Ol bond 
lengths obtained from a refinement with isotropic param- 
eters are reported in Table |TJ Apparently there is no OIE 
on atomic positions within the resolution of the diffrac- 
tometer and the fitting procedure. The results obtained 
using an anisotropic refinement are identical within the 
range of error. 

In the following we express the OIE on a quantity x 
such as lattice parameters a, b, c etc. by Ax = x 7if — 
(a; 16 and x 18 are short for the quantity x measured 
in La2Cu 16 04 and La2Cu 18 04, respectively). The re- 
sults for OIE on the lattice parameters and the Cu-Ol 
bondlength at 15 and 290 K are summarized in Table HT1 
the OIE turned out to be in the range of 0.01 % for 
all the parameters. Interestingly the orthorhombicity 
is more pronounced in the 18 O sample: we found that 
Oi 6 /&i6 = 0.98435(1) at T = 15 K and 0.99101(1) at 
T = 290 K whereas a 18 /b 18 = 0.98419(1) at 15 K and and 
0.99093(1) at 290 K. The Cu-Ol bondlength vs. temper- 
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FIG. 3: Cu-Ol bond lengths vs. temperature T for 
La2Cu 16 C>4 (solid circles) and La2Cu 18 C>4 (open circles). 



ature is displayed in Fig. [3] 

In Fig. |4] the temperature evolution of the isotropic 
parameters Bi SO for Cu and 01 is illustrated for both 
the 16 O and 18 O sample. The B iso parameters increases 
with temperature, but below temperatures of approxi- 
mately 100 K a saturation sets in. We note that the 
points at 15 K and at 290 K are measured with high 
statistics as opposed to the intermediate points which 
explains the larger error bars and the offset when going 
from the high-statistics to the low-statistics data. The re- 
sults from the high-statistics data are summarized in Ta- 
ble lllll The magnitude of the J3i SO -parameters scale qual- 
itatively with the inverse mass of the atoms, more pre- 
cisely B lso (02) > B iso (01) > B iso (Cu) > B lso (La) at 
15 K and B lso (02) > B lso {01) > B lso (La) > B lso {Cu) 
at 290 K, respectively. 

The corresponding OIE on _Bi SO -parameters are given in 
Table HV1 We found that Bi SO of the oxygen atoms 01 
and 02 exhibit a negative OIE, whereas the OIE is pos- 
itive for La and Cu. 

The results of the anisotropic thermal parameters are 
also reported in Table IHT1 Since data were collected from 
a powder sample it is only possible to refine the diagonal 
elements Un, U22 and C/33. From the order of magni- 
tude of Ua, 1 < i < 3 we deduce in which direction the 
movements are the largest: at 15 K U max = U22 for La, 
Cu and 01 and U max = U n for 02. At 290 K we find 
U max = U22 for La, U max = Un for 02 and U rnax = U33 
for Cu and 01. Similarly to the isotropic refinement the 
oxygen atoms showed a negative OIE, whereas La and 
Cu atoms had a positive OIE in all the directions. In 
particular we found that the OIE on Un is huge for La 
(~ 90 %)and Cu (~ 47 %) at 15 K. 




FIG. 4: Temperature dependence of the isotropic thermal 
parameters Bi ao for Cu and Ol. 



IV. DISCUSSION 

Our neutron diffraction study showed that the samples 
La 2 Cu 16 04 and La 2 Cu 18 04 crystallized both in the or- 
thorhombic space group Bmba, which is well known in 
literature^ 7 .. We would like to emphasize that they are 
single-phased which confirms our results of the analysis 
of the oxygen content by means of hydrogen reduction 
which yielded an oxygen content of 4.004 ± 0.005. The 
small excess oxygen concentration explains the low val- 
ues of the peak in the susceptiblity which is related to the 
Neel temperature Tjy*. The oxygen isotope substitution 
raises the peak in the susceptibility by 2.4 ± 0.2 K as op- 
posed to results reported in RefAi. However, we would 
like to point out that Zhao et al^ found that the values 
of TV and also the OIE on depends strongly on the 
preparation and annealing conditions of the samples. 
The structural parameters obtained from our neutron 
diffraction measurements on La2Cu 16 04 are in excellent 
agreement with literatur o 17 ' 18 . Measurements on HRPT 
which is very well suited to determine lattice constants 
with highest accuracy revealed a negative OIE on the 
lattice constants, see Table HU We found a tendency for 
the Cu-Ol bond lengths to be larger in the 16 sample 
as compared to the 18 sample by the order of 10 -4 , see 
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a (A) b(A) c (A) 


y(La) «(La) z(01) y{02) z(02) 


Cu-Ol (A) 


15 K 

16 
18 G 


5.33305(2) 5.41783(3) 13.10260(7) 
5.33185(3) 5.41751(3) 13.09979(8) 


-0.00839(13) 0.36149(4) -0.00832(6) 0.04101(12) 0.18302(8) 
-0.00855(14) 0.36158(4) -0.00834(6) 0.04110(13) 0.18293(8) 


1.90369(5) 
1.90343(5) 


290 K 

16 
i8 G 


5.35479(3) 5.40338(3) 13.14810(9) 
5.35385(3) 5.40283(3) 13.14614(9) 


-0.00683(19) 0.36135(4) -0.00746(8) 0.03479(17) 0.18308(9) 
-0.00655(20) 0.36138(4) -0.00736(8) 0.03463(17) 0.18291(9) 


1.90434(5) 
1.90401(5) 



TABLE I: Lattice parameters, atomic coordinates and the Cu-Ol bond lengths for La2Cu O4 and La2Cu O4 at 15 and 
290 K as obtained in an isotropic refinement. In the nonstandard space group Bmba (isomorphic to the standard Cmca, 
No. 64) used in the present work, the atomic positions are La [0, y, z], Cu [0,0,0], Ol [1/4, 1/4,2], 02 [0, y,z]. 



T 


Aa [10~ 4 ] 


Ab [10- 4 ] 


Ac [10~ 4 ] 


A Cu-Ol [10~ 4 ] 


15 K 


-2.25(9) 


-0.59(11) 


-2.14(14) 


-1.37(53) 


290 K 


-1.76(13) 


-1.02(11) 


-1.37(53) 


-1.73(53) 



TABLE II: OIE on the lattice parameters in La2Cu04 as ob- 
tained by an isotropic refinement. 
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FIG. 5: Lattice parameter a as a function of temperature 
with displacement due to thermal motion as deduced from 
Biso parameters. 



Table HI1 No OIE could be detected regarding the atomic 
positions within the instrumental resolution and the ac- 
curacy of the Rietvield refinement. 

The values of the isotropic and anisotropic thermal pa- 
rameters Bi SO and Uu, 1 < i < 3 are in agreement with 




FIG. 6: Anisotropic thermal parameter for Cu in the 1 O 
(gray) and 18 O (white) sample at 15 K (left panel) and 290 K 
(right panel). 



results reported in Refs . 17 ' 18 and in Rcf. 19 , respectively. 
We found a significant OIE on the zero point motion 
in terms of the thermal parameters up to approximately 
50 % at 15 K. The OIE still exists at room temperature, 
however it is substantially lower, its maximum is ~ 10 %, 
see Table IIV1 Interestingly, we observed a negative OIE 
on the thermal parameters - both B iso and Uu , 1 < i < 3 
- for Ol and 02, whereas the OIE is positive in case of 
La and Cu which is qualitatively consistent with the sum 
rule for lattice vibrations^. Furthermore we would like 
to emphasize that the sign of the OIE on the anisotropic 
thermal parameters Uu is identical for the i = 1,2,3 di- 
rection with the sign of the OIE on the Bi SO parameters 
for each atom. 

Interestingly the OIE on the anisotropic thermal parame- 
ters varies significantly for the different direction of ther- 
mal motion at 15 K: AUu « 0.9 down to AU22 ~ 0.3 
and AUn ~ 0.5 down to AU33 « 0.2, for La and Cu 
respectively. In the case of the oxygen atoms Ol and 
02 the order of magnitude of the OIE on the anisotropic 
thermal parameter stays in a comparable range ~ 0.2. 
At 290 K the changes in the OIE are substantial for La 
and Cu: while the OIE on Uu (E/33) vanishes within the 
present accuracy, we obtained AU33 « 0.2 (AUu ~ 0.3) 
for La (Cu). The direction where the OIE is most pro- 
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15 K 


290 K 




Biso U\\ U22 


Bi S o U\\ U22 £^33 


10 La 

18 o 


0.075(6) 0.00059(16) 0.00171(18) 0.00060(20) 
0.115(6) 0.00113(18) 0.00218(19) 0.00111(32) 


0.367(7) 0.00526(17) 0.00553(21) 0.00371(18) 
0.407(7) 0.00540(19) 0.00612(22) 0.00458(18) 


10 Cu 

18 o 


0.109(8) -0.00019(20) 0.00214(23) 0.00212(35) 
0.144(9) 0.00028(28) 0.00263(30) 0.00248(35) 


0.320(10) 0.00145(26) 0.00386(28) 0.00764(35) 
0.360(10) 0.00192(18) 0.00443(30) 0.00808(35) 


lb 01 

18 o 


0.234(9) 0.00237(26) 0.00330(25) 0.00289(35) 
0.183(9) 0.00176(27) 0.00256(27) 0.00235(35) 


0.520(11) 0.00374(25) 0.00480(27) 0.01207(36) 
0.471(11) 0.00357(26) 0.00430(28) 0.01083(36) 


i0 02 

18 o 


0.416(9) 0.00693(24) 0.00431(31) 0.00389(26) 
0.339(10) 0.00601(24) 0.00324(31) 0.00302(26) 


0.970(13) 0.01923(35) 0.01276(44) 0.00498(26) 
0.915(13) 0.01862(36) 0.01118(44) 0.00456(35) 



TABLE III: B iso -parameters in A 2 and [/^-parameters (1 < i < 3) in A 2 of La 2 Cu0 16 4 and La 2 Cu 18 04 at 15 and 290 K for 
an isotropic and anisotropic refinement, respectively. 





15 K 


290 K 




ABi SO 


ABi SO 


La 


0.53(16) 


0.11(4) 


Cu 


0.32(15) 


0.13(6) 


01 


-0.22(8) 


-0.094(4) 


02 


-0.19(4) 


-0.06(3) 



TABLE IV: OIE on isotropic thermal parameters ABi ao for 
La 2 Cu0 4 at 15 and 290 K. 

nounced is not identical with the direction where Uu is 
largest. 

In order to visualize the effect of the isotropic thermal 
parameters we show in Figure [5] the temperature evolu- 
tion of the lattice parameter together with the mean dis- 
placement a±ao(Cu) deduced from the isotropic thermal 
parameters for the copper atom Bi SO {Cu) 

a ° = {^%- (1) 

In this scale it is striking that the OIE on the lattice pa- 
rameter a is not visible anymore, even the temperature 
dependence of a ± a as a whole expression is small. But 
the OIE ona±ao due to the movement of the copper ion 
can be observed in terms of a being subject to a larger 
variation due to the thermal motion of Cu in the 18 O 
compound than in the 16 compound. 
Figure [5] illustrates the anisotropic thermal parameters 
for Cu in the shape of an ellipsoid with axes Uu, U22 
and C/33 corresponding to thermal motion. At 15 K the 
zero-point motion is almost isotropic in the y-z direction, 
while there is a considerable anisotropic deformation at 
room temperature. For 01 which is also in-plane the 
anisotropy is similar, the elongation in the z-direction 
is also most pronounced at 290 K while the zero point 
motion is almost isotropic. The out-of plane ions 02 
and La behave in the different way: at 15 K La is al- 
most isotropic in the x-z direction and 02 is almost 
isotropic in the y-z direction, whereas the maximum dis- 
placement is along the y- and z-direction for La and 02, 
respectively. At room temperature the thermal motion 



is highly anisotropic and the largest elongation is along 
the x-direction for both out-of-plane ions La and 02. 



V. CONCLUSIONS 

In the present work we have studied the influence 
of the oxygen isotope effect ( 18 vs. 16 78%) 
on structural and thermal parameters by means of 
high-resolution neutron diffraction experiments on the 
antiferromagnetic insulator La2Cu04. The oxygen con- 
centration in the 16 O and 18 O samples was determined 
to be 4.004 ± 0.005 by hydrogen reduction. The excess 
oxygen content was confirmed to be small by neutron 
diffraction measurements which revealed the samples 
to be single-phased (Cmca). We found a nonvanishing 
negative OIE on the lattice parameters in the order of 
0.01%. The OIE on the Cu-Ol bond length turned out 
to be negative as well. We found a considerable OIE 
on the zero point motion in terms of the isotropic and 
anisotropic thermal parameters (up to 100% for La), it is 
positive for the atoms La and Cu, whereas it is negative 
for 01 and 02. A refinement using anisotropic thermal 
parameters yielded qualitatively the same result, in 
particular the sign of the OIE on the anisotropic thermal 
parameters is identical for the x-, y- and z-direction 
with the sign of the OIE on the isotropic thermal 
parameters. Moreover, we discovered that the OIE on 
the anisotropic thermal parameters varies significantly 
with the direction of thermal motion e.g. AUn « 0.9 
and MJn w 0.3 for La at 15 K. The huge OIE on the 
isotropic thermal parameters can therefore be attributed 
mainly to the thermal motion in the x-direction. 
Furthermore our diffraction study revealed that the 
anisotropic deformation of the ellipsoid of thermal 
motion is different for in-plane and out-of-plane ions. 
The elongation is most pronounced along the z-direction 
for Cu and 01 (in-plane) and along the x-direction for 
La and 02 at 290 K. Zero point motion is isotropic 
at least in two directions. The preferred direction of 
thermal motion is identical for the 16 and 18 sample. 
Thermal motion in La2Cu0 4 is in the range of 10-20 
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meV 21 whereas the energy scale of the antiferromagnetic 
exchange J is of the order of 100 meV— . Thus the 
displacement of the ions due to thermal motion is 
almost static as compared to the electronic scale of the 
hopping integral t and the antiferromagnetic exchange J. 

We hope that the present detailed structural data 
will stimulate theoretical studies on how the OIE on 
structural parameters as well as on the zero point motion 
affects the hopping integral t and the antiferromagnetic 



exchange J in La2Cu04. 
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